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摘要 : 高 空 云层 导致 所 观测 的 Ha 全 日 面 像 上 履 盖 有 一 层 云 污染 ， 使 得 图 像 上 的 太阳 活 
动 细节 变 得 模糊 不 清 。 为 了 能 够 实时 探测 云 污染 ， 并 及 时 显示 修复 后 的 图 像 ， 采 用 图 形 处 理 
器 技术 实现 了 一 个 Ha 全 日 面 云 污染 实时 识别 和 修复 系统 。 该 系统 主要 在 统一 计算 设备 架构 
(Compute Unified Device Architecture, CUDA) 环境 下 利用 图 形 处 理 器 并 行 实现 : (1) 二 值 化 图 
像 椭圆 长 短 轴 比 值 法 识别 重度 云 污 染 图 像 ， (2) 临 边 绰 瞳 曲 线 中 心 对 称 法 识别 可 修复 云 污 染 
Bik; (3) 频 域 巴特 沃 斯 低 通 滤波 法 去 除 云 污染 。 通 过 对 系统 中 各 运算 在 图 形 处 理 器 中 花费 
的 时 间 进 行 详细 测量 ， 发 现 傅 里 叶 正 反 变 换 和 频 域 滤波 占用 了 图 形 处 理 器 总 处 理 时 间 的 
52. 9% ， 是 系统 中 最 耗 时 的 。 然 而 ， 相 对 于 1 min 的 观测 时 间 间 隔 ， 约 0.7s 的 总 处 理 时 间 可 
以 满足 实时 显示 的 需要 。 另 外 ， 通 过 对 修复 后 的 图 像 做 质量 评价 ， 验 证 了 所 采用 的 修复 算法 
可 以 有 效 地 去 除 云 污染 ， 并 对 太阳 活动 细节 影响 较 小 。 最 后 讨论 了 系统 存在 的 不 足 和 需要 进 
行 的 改进 。 
关键 词 : 云 污染 ， 并 行 计算 ; 识别 和 修复 ; 巴特 沃 斯 低 通 滤波 器 
中 图 分 类 号 : P111; TP274 ”文献 标识 码 : A ”文章 编号 ; 1672-7673(2016)02-0242-08 


Ha 全 日 面 观测 为 监测 太阳 色 球 层 活动 现 象 提供 大 量 的 数据 。 全 球 拥有 众多 的 Ha 全 日 面 观测 
站 ， 如 全 球 日 震 网 ( Global Oscillation Network Group, GONG) HA EER 6 个 观测 点 中 ,国内 的 北京 怀 
柔 观 测 基地 等 等 。 虽 然 各 站 点 在 选 址 时 进行 了 严格 的 筛选 ， 但 是 实际 观测 中 会 出 现 大 量 的 多 云天 
气 呈 。 在 这 样 的 天 气 观 测 ， 经 常 出 现 太 阳 局 部 被 云 履 盖 的 情况 ， 导 致 图 像 细 节 模 糊 。 图 1 显示 了 
Cerro Tololo 美洲 际 天 文 台 在 2015 年 4 月 20 日 观测 的 3 张 图 像 ， 可 以 看 出 后 两 张 图 像 存 在 不 同 程度 的 
云 污 染 现 象 。 

云 覆 盖 是 导致 全 日 面 像 畸 变 的 一 种 情况 。 为 了 检测 和 校正 各 种 畸变 ， 前 人 提出 了 一 些 关 于 Ha 全 
日 面 像 的 质量 评价 和 修复 算法 。 如 文 [3] 采 用 双 尺 度 滤波 归 一 化 图 像 强 度 ; 文 [4] 根 据 图 像 临 边 昏 暗 
曲线 的 峰 度 判断 图 像 质 量 ; 文 [5] 根 据 最 小 相关 系数 法 判断 图 像 是 否 有 云 污染 ; 文 [6-7] 采 用 大 尺度 
中 值 滤波 修复 畸变 图 像 和 去 除 云 污染 , 文 [5] 使 用 多 尺度 形态 学 滤波 法 修复 云 污染 图 像 等 。 这 些 方法 
主要 为 了 在 后 续 的 太阳 物理 研究 中 ， 能 在 观测 图 像 上 更 好 地 切割 、 识 别 和 测量 各 种 太阳 活动 特征 结 
构 ， 如 炊 斑 、 暗 条 等 等 。 因 此 ， 并 不 涉及 执行 效率 和 运行 速度 的 问题 。 

如 果 能 将 这 些 方法 用 于 实时 探测 是 否 有 云 污染 ， 并 实时 显示 修复 后 的 图 像 ， 将 非常 有 利于 观测 的 
进行 和 观测 质量 的 判断 。 但 这 必然 要 求 相应 的 优化 算法 和 计算 机 的 快速 处 理 能 力 。 为 此 基于 文 [5] 的 
处 理 方法 ， 采 用 图 形 处 理 器 技术 实现 了 一 个 实时 全 日 面 像 云 污染 的 识别 和 修复 系统 。 在 硬件 环境 下 ， 
处 理 一 张 云 污染 图 像 消耗 的 总 时 间 约 为 0.7 s。 根 据 图 形 处 理 吉 技术 的 特点 ， 采 用 频 域 巴特 沃 斯 滤波 
器 ， 并 使 用 图 像 相 似 度 算法 对 无 云图 像 和 云 污染 修复 后 的 图 像 进 行 局 部 相似 度 检 验 。 结 果 显 示 ， 在 满 
足 处 理 速 度 的 同时 ， 这 一 滤波 算法 可 以 在 很 大 程度 上 保留 太阳 的 活动 细节 ， 满 足 观测 实时 显示 的 需要 。 
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(a) 
Al (a), (b), (c) 3 MARA Cerro Tololo 美洲 际 天 文 台 分 别 于 当地 时 间 
2015 年 4 月 20 日 (16:07、16:08、20:17) 拍摄 
Fig.1 (a)，(b) and (c) are three Ha full-disk solar images, which were observed 


by the GONG Ha instrument of the CTIO on April 20, 2015 


1 算法 原理 分 析 


(c) 


对 云 污染 图 像 的 处 理 包括 两 步 : 云 污染 的 检测 和 云 污染 的 修复 。 其 中 检测 包括 : (1) 判 断 全 日 面 
像 是 否 为 一 个 圆 ， 否 则 认为 是 重度 云 污染 的 图 像 而 放弃 修复 ; (2) 判断 全 日 面 像 的 临 边 借 暗 曲线 是 否 
圆 对 称 ， 如 果 不 对 称 ， 说 明 有 云 污染 。 云 污染 的 修复 采用 文 [5] 提 出 的 模板 法 ， 即 通过 和 一 个 标准 全 
日 面 像 模板 进行 对 比 ， 并 通过 滤波 得 到 云层 透 过 率 ， 从 而 修复 被 污染 的 图 像 。 


1.1 重度 云 污染 图 像 的 判断 

正常 情况 下 的 全 日 面 像 其 太阳 圆 面 与 天 空 青 
景 的 强度 值 相 差 很 大 ， 可 以 用 一 个 高 于 天 空 背景 
的 国 值 对 其 进行 二 值 化 。 当 云 污染 过 重 时 ， 日 面 
像 就 会 残缺 。 图 2 为 图 1(c) 二 值 化 后 的 图 像 。 

因此 ， 可 先 将 二 值 化 图 像 假 定 为 一 个 椭圆 ， 
通过 计算 其 长 轴 与 短 轴 的 比值 环 判 断 是 否 残 缺 。 
对 良好 观测 条 件 下 的 GONG 图 像 统计 表明 ， 其 五 
值 为 1.0010 + 0.0015。 而 图 2 W E (EN 1.938, 
可 判定 为 重度 污染 的 图 像 。 在 实际 系统 中 ， 凡 五 
值 高 于 1.1 的 都 归 为 重度 污染 的 图 像 。 
12 无 云 污染 图 像 的 判断 

对 于 二 值 化 后 为 一 个 圆 的 图 像 ， 多 数 没有 云 
污染 。 根 据 文 [5] 的 方法 ， 可 通过 判断 临 边 昏 暗 
曲线 是 否 中 心 对 称 判断 是 否 有 云 污染 存在 。 主 要 
步骤 包括 : (1) 将 太阳 图 像 中 心 化 并 将 其 从 直角 
坐标 系 转 到 极 坐 标 系 ; (2) 获 取 极 坐标 图 像 等 分 


图 2 图 1(c) 二 值 化 处 理 后 图 像 
Fig.2 Binary image processing result of Fig. 1(c) 


4 部 分 中 值 临 边 昏 暗 曲线 ; (3)4 AHL SF HH A LOR RL, ECB Bc AR RL, 


如 果 无 云 污染 ， 这 个 系数 接近 1， 反映 了 临 边 昏暗 


线 的 圆 对 称 性 。 当 有 云 污染 时 ， 这 个 系数 会 


变 小 。 在 系统 中 采用 的 阔 值 为 0.95。 图 3 显示 了 图 1(b) 的 4 条 临 边 昏 暗 曲线 ， 其 最 小 相关 系数 为 
0. 32， 而 相应 的 图 1(a) 为 0.99。 据 此 判定 图 1 中 (a) 不 存在 云 污染 ，(b) 存 在 可 修复 的 云 污 染 。 


1.3 云 污 染 的 修复 


云 污 染 全 日 面 像 的 成 像 模型 可 简化 为 以 下 公式 : 
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T(x, y)=c(x, y)t(x, y), 


EH, c(a, y) 为 云层 透射 率 ; (a, y) 为 没有 云 Regen 

污染 的 全 日 面 像 ，7(x, y) 为 实际 观测 的 云 污染 ple Sy, 
全 日 面 像 。 修 复原 理 可 简单 表示 成 1(%, y)= T(x， a) 
y)/e(x, y) 。 修 复 的 步骤 为 ， (1) 选 择 一 张 无 去 图 8 5 \ 
像 (如 图 1(a) AEREE RRR T i 

4(a) ) 作 为 模板 ; (2) 用 有 云 污染 的 图 像 (图 1 the second quadrant 

Cb) ) Ban FE AE COA, FSA KH oo 

动 细节 的 云层 透射 率 (图 4(b)); (3) 将 图 4(b) o F i 
转换 到 频 域 并 利用 巴特 沃 斯 滤波 器 对 其 进行 低 通 r 

滤波 ， 扣 除 太阳 活动 细节 的 影响 ,得 到 云层 透射 图 3 图 1(b) 对 应 的 四 个 象限 的 中 值 临 边 措 上 暗 曲 线 
率 ( 图 4(c)); (4) 将 原始 的 云 污染 图 像 除 以 云层 Fig.3 The median limb-darkening profiles of four 
透射 率 得 到 最 终 去 除 云 污染 的 图 像 (图 4(d) ) 。 eee eae Ae) 


(a) (b) (c) (d) 
图 4 (a) HA l(a) PARE RE; (p) 含 太阳 活动 细节 的 云层 透射 率 ; 
(ce) 去 除 太阳 活动 细节 的 云层 透射 率 ; (d) 去 除 云 污染 的 图 像 
Fig.4 (a) The standard profile of Fig. 1(a); (b) The clouds transmittance with solar features; 


(c) The clouds transmittance without solar features; (d) The solar image without cloud pollution 


在 第 (3) 步 发 现 ,图 4(b) 中 的 太阳 活动 细节 在 频 域 图 中 主要 分 布 于 以 中 心 为 原点 ， 半 径 为 D, 的 
圆 外 区 域 。 因 此 在 巴特 沃 斯 滤波 器 中 ， 只 需 选 择 合适 的 半径 Du( 截 止 频 率 ) 将 其 圆 外 的 高 频 部 分 滤 除 ， 
即 可 在 滤波 后 得 到 扣除 太阳 活动 细节 的 云层 透射 率 。 经 过 统计 不 同时 段 的 云 污染 图 像 ， 将 系统 中 截止 
频率 取 为 图 像 尺 寸 的 0. 003 倍 。 

为 充分 发 挥 图 形 处 理 器 的 并 行 处 理 优势 ， 采 用 频 域 滤波 法 。 为 此 ， 采 用 图 像 结构 相似 度 ( Structural 
Similarity Index Measurement, SSIM) 算法 证 实 修复 的 可 靠 性 。 选 两 幅 拍 摄 时 间 仪 隔 一 分 钟 的 图 像 ， 
第 1 幅 是 干净 的 图 像 ， 第 2 幅 为 云 污染 的 图 像 ， 对 第 2 幅 图 像 进 行 修 复 得 到 第 3 幅 图 像 。 图 5 显示 了 
这 3 幅 图 像 的 局 部 区 域 。 计 算得 到 无 云图 像 ( 图 5(a)) 与 有 云图 像 (图 5(b)) 的 图 像 结 构 相 似 度 为 
0. 2997， 无 云图 像 (图 5(a) ) 与 去 除 云图 像 ( 图 5(c)) 的 图 像 结构 相 似 度 为 0.755 8。 可 以 看 出 去 除 云 
污染 之 后 图 像 的 结构 相似 性 明显 高 于 去 除 云 污 染 之 前 的 图 像 ， 说 明 频 域 滤 波 也 可 以 相当 程度 地 去 除 云 
污染 ， 并 保留 了 太阳 的 活动 细节 。 


2 云 污染 图 像 识 别 和 恢复 系统 的 图 形 处 理 带 实现 


形 处 理 器 的 设计 初 囊 是 为 了 将 计算 机 中 关于 图 形 图 像 的 处 理 和 显示 从 中 央 人 处理 器 中 独立 出 来 。 
由 于 图 形 处 理 器 采用 了 单 指令 多 数据 流 ' 的 数据 操作 方式 ， 因 而 并 行 计 算 能 力 很 强 。 为 了 将 图 形 处 
理 器 的 多 线程 并 行 能 力 运 用 到 通用 计算 上 来 ，NVIDIA 公司 推出 了 统一 计算 设备 架构 这 一 编程 模 
HY) 。 全 日 面 像 云 检测 和 修复 算法 具备 数据 计算 密度 高 、 可 并 行 性 强 的 条 件 ， 因 而 十 分 适合 用 图 形 


X N 


201711.01098v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


2 期 张 能 维 等 ，Ha 全 日 面 云 污染 实时 识别 和 修复 系统 245 


人 处理 器 处 理 。 并 且 CUDA 中 有 很 多 现成 的 函数 库 ， 如 CUBLAS 库 中 的 大 量 针 对 和 矩阵 的 基本 运算 函数 
以 及 CUFFT 库 中 的 傅 里 叶 正 反 变 换 函 数 ， 可 以 极 大 提高 程序 的 运行 效率 。 例 如 ,， 文 [11] 利用 图 形 处 
理 器 技术 在 国家 天 文 台 怀 柔 观 测 站 实现 一 个 实时 深 积分 磁场 观测 系统 ,， 文 [12] 利 用 图 形 处 理 器 技术 
实现 了 抚 仙 湖 1 m 真空 太阳 望远镜 level 1 级 图 像 选 帧 。 


=r 


(a) (b) O) 
图 5 (a), (b), (ce) 分 别 为 图 1(a)、 图 1(b)、 图 3(d) 自 中心 截 取 的 1200x1200 像 素 局 部 图 
Fig.5 (a), (b) and (c) are three 1200 x 1200 partial image, which are intercepted 
from the center of Fig. 1(a), Fig. 1(b) and Fig. 3(d) 


图 6 展示 了 系统 流程 图 。 该 系统 首先 要 显示 实时 处 理 界面 ， 然 后 循环 检测 是 否 有 更 新 的 图 像 。 如 
果 发 现 更 新 图 像 先 显示 原 图 ， 再 做 云图 像 检 测 和 修复 处 理 ， 最 后 把 处 理 后 的 图 像 显 示 在 界面 上 与 原 图 
像 对 比 。 如 果 没 有 检测 到 更 新 图 像 ， 则 界面 一 直 显 示 之 前 图 像 的 处 理 结果 。 

其 中 中 央 处 理 器 所 要 完成 的 主要 任务 包括 : (1) 将 原始 FITS 格式 图 像 数 据 从 硬盘 读 和 中央 处 理 
器 内 存 并 显示 在 界面 上 ; (2) 将 图 像 数据 从 中 央 处 理 需 内 存 复 制 到 图 形 处 理 器 显存 ; (3 ) 将 处 理 后 得 
到 的 数据 从 图 形 处 理 器 显存 复制 回 中 央 处 理 器 内 存 并 显示 。 

图 形 处 理 需 执行 系统 的 核心 任务 一 一 云 检 测 和 去 去除。 将 图 形 处 理 器 中 计算 任务 大 致 分 为 两 类 
一 类 是 通过 直接 调用 CUDA 库 函 数 实现 。 另 一 类 是 自行 编写 的 利用 线程 并 行 实现 的 函数 ， 其 执行 的 
主要 思想 为 : (1) 根 据 数 和 矩阵 点 的 总 数 分 配 总 线程 数 ， 通 过 线程 号 索引 让 和 矩 阵 中 每 个 点 与 线程 一 一 对 
应 ; (2) 根 据 判断 条 件 将 每 个 线程 要 完成 的 计算 任务 分 类 ， 同 一 判断 条 件 下 所 有 线程 的 计算 任务 相 
同 ; (3) 所 有 线程 并 行 执行 各 自分 配 的 计算 任务 。 

在 检测 重度 云 污染 图 像 环 节 ， 首 先 并 行 实现 图 像 二 值 化 处 理 。 椭 圆 长 轴 与 短 轴 并 行 求 取 方法 为 ; 
(1) 并 行 求 取 椭 圆 重心 位 置 ，(2) 并 行 得 到 椭圆 上 各 点 与 其 重心 的 协 方差 。 最 后 将 所 得 协 方差 元 素 代 
入 公式 即 可 得 到 椭圆 长 短 轴 之 比 。 据 此 如 果 判 断 图 像 云 污染 过 重 ， 则 丢弃 ， 否 则 继续 下 面 步骤 。 

接 下 来 要 检测 图 像 是 否 存 在 可 修复 的 云 污染 。 首 先 并 行 实现 将 图 像 上 的 全 日 面 从 日 心 (对 应 上 面 
所 求 重心 ) 位 置 平移 到 图 像 中 心 位 置 ， 即 得 到 中 心 化 图 像 。 中 心 化 过 后 ， 再 并 行 实现 用 双 线 性 插值 法 
将 图 像 转换 至 极 坐标 系 下 。 紧 接着 对 极 坐标 图 像 角度 的 4 个 等 间距 区 域 分 别 进行 排序 ， 排 序 方法 是 并 
行 实现 的 双 调 排序 算法 。 从 排 完 序 的 极 坐标 图 像 中 直接 取出 所 需 的 4 条 临 边界 上 暗 曲 线 ， 将 两 条 曲线 用 
现 有 函数 进行 点 积 、 求 和 等 运算 ,再 将 结果 代入 公式 即 可 得 两 条 曲线 的 相关 系数 。 利 用 比较 所 得 最 小 
相关 系数 如 判断 有 云 污染 ， 则 不 需要 修复 ， 否 则 继续 下 面 去 除 云 污染 的 步骤 。 

在 使 用 模板 法 去 除 云 污染 时 ， 要 用 到 模板 的 标准 临 边 昏 暗 轮廓 。 生 成 模板 算法 与 检测 云 算 法 很 相 
似 ， 图 7(a) 展 示 了 Ho 全 日 面 像 生成 模板 系统 的 界面 。 先 通过 并 行 相 除 得 到 带 有 太阳 活动 细节 的 云 
层 透 射 率 。 紧 接着 利用 CUFFT 库 中 函数 进行 傅 里 叶 变 换 将 图 像 转 换 到 频 域 。 再 并 行 生成 一 个 巴特 沃 
斯 低 通 滤波 器 ， 将 其 与 频 域 图 像 并 行 相 乘 ， 即 可 得 到 滤 除 高 频 部 分 的 频 域 图 像 。 然 后 将 滤 除 高 频 部 分 
的 频 域 图 像 利用 传 里 叶 逆 变 换 转 回 时 域 得 到 消除 太阳 活动 细节 的 云层 透射 率 。 最 后 将 中 心 化 云 污染 像 
与 云层 透射 率 并 行 相 除 得 到 最 终 去 除 云 污染 的 图 像 。 

对 于 系统 中 一 些 重要 参数 ， 根 据 以 往 的 数据 实验 统计 设置 默认 值 。 为 了 让 用 户 可 以 在 合理 范围 内 进 
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行 调整 ， 同 时 将 这 些 参数 框 放 在 界面 上 。 如 图 7(b) 为 Ha 全 日 面 像 云 污 染 实时 识别 和 修复 系统 界面 。 


图 像 数 据 二 值 化 


根据 二 值 化 医 
椭圆 长 短 轴 2 


hl ae 
新 图 像 传 


显示 前 一 由 
图 像 处 理 结果 PES 
FERRI Ea 

平移 到 图 像 | 


读 入 图 像 数 据 并 
显示 在 界面 上 


将 极 坐 标 图 像 排序 并 
求 其 最 小 相关 系数 


图 像 
作为 最 终结 果 


图 6 Ha 全 日 面 云 污染 实时 识别 和 修复 系统 
Fig.6 The flow chart of the real-time processing system for detecting and 


removing cloud shadows on Ha full-disk solar image 
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(a) (b) 
图 7 (a) 生 成 模板 系统 界面 ;(b) 云 污染 实时 识别 和 修复 系统 界面 
Fig.7 (a) The system interface of generating templates; (b) The system interface 


of real-time detecting and removing cloud shadows 
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3 系统 性 能 测试 


实验 采用 的 计算 机 配置 为 : 中 央 处 理 器 为 英特尔 Pentium( 奔腾 ) 双核 E5200 处 理 器 ， 图 形 处 理 器 
卡 为 Nvidia 公司 的 GeForce GT 430, 操作 系 统 为 Windows 7 旗舰 版 64 位 。 软 件 开发 环境 为 Microsoft 
Visual Studio 2010 下 基于 对 话 框 的 微软 基础 类 库 ( Microsoft Foundation Classes, MFC) 和 CUDA 6.5, 

对 系统 中 各 主要 运算 在 图 形 处 理 器 中 花费 的 时 间 进 行 测量 ， 测 试 数据 为 图 1(b) 云 污染 图 像 ( 尺 十 
为 2 048x2 048 (RA), ， 得 到 如 表 1 的 测试 结 
表 1 系统 中 各 运算 在 图 形 处 理 器 中 处 理 时 间 


Table 1 The processing time of each operation in GPU 
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运算 任务 运算 时 间 /ms 占 总 时 间 比 /% 
将 图 像 原始 数据 从 中 央 处 理 器 复制 到 图 形 处 理 器 14. 2 3.8 
图 像 二 值 化 3.9 
重度 云 污染 求 各 像素 点 横 纵 坐标 5.8 
图 像 的 判断 求 椭圆 重心 坐标 4.1 33.2 8.8 
求 椭 圆 各 点 与 重心 偏 移 量 6.9 
求 椭圆 长 短 轴 比 值 4.8 
图 像 中 心 化 10.6 
于 图 像 直角 坐标 到 极 坐 标 转换 8.4 
无 云 污染 eas 
图 像 的 判断 极 坐 标 图 像 行 排序 26. 8 48. 8 13.0 
Puh 4 AFP LM ido es HA 0.3 
计算 4 条 曲线 的 最 小 相关 系数 1.4 
云 污 染 图 像 除 以 模板 标准 临 边 昏 上 暗 轮 局 9.9 
将 实数 转 为 复数 2.5 
傅 里 叶 变 换 77.5 
云 污 染 图 像 生成 巴特 沃 斯 低 通 滤波 器 45. 6 
的 修复 低 通 滤 波 3.5 
逆 傅 里 叶 变换 72.3 
复数 转 为 实数 9.4 
云 污染 图 像 除 以 云层 透射 率 9.9 
将 处 理 结 果 数据 从 图 形 处 理 器 复制 回 中 央 处 理 器 21.7 5.7 
其 他 17. 5 4.7 
图 形 处 理 需 处 理 总 时 间 376. 1 100 


经 测试 , 除了 图 形 处 理 器 的 处 理 时 间 外 ， 中 央 处 理 器 的 总 运行 时 间 为 342 mso HF, PRAHA 
花费 的 时 间 包 括 : 从 硬盘 读 取 原始 图 像 数 据 和 模板 图 像 数 据 (154 ms) 、 在 界面 上 显示 原始 图 像 和 修复 
后 的 图 像 (160 ms) 、 其 余 运 算 (28 ms) 。 所 以 ， 该 系统 处 理 总 时 间 为 718 ms, 


4 讨 论 


图 形 处 理 器 对 高 密度 图 像 数据 的 并 行 处 理 能 力 在 系统 中 得 到 了 充分 体现 。 对 于 一 幅 2 048x2 048 


像素 的 云 污染 Ha 全 日 面 图 像 ， 系 统 在 图 形 


处 理 融 中 完成 从 检测 到 修复 等 处 理 运 算 的 总 时 间 为 376. 1 mso 


男 外 ， 系 统 在 中 央 处 理 屁 中 完成 读 图 像 、 显 示 图 像 等 处 理 总 耗 时 为 342 ms。 针 对 常规 1 min 采样 间隔 


的 Ha 全 日 面 观测 ， 可 以 满足 观测 过 程 中 对 处 理 结果 实时 


a 


BY 
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势 ， 在 修复 过 程 中 采用 频 域 巴 特 沃 斯 低 通 滤波 的 方法 。 通 过 使 用 图 像 结构 相 似 度 算法 对 修复 后 图 像 做 
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质量 评价 ， 发 现 该 方法 与 前 人 所 用 的 中 值 滤波 法 和 多 尺度 形态 学 滤波 法 修复 效果 相当 。 系 统 除 了 可 以 
处 理 GONG 的 数据 ， 还 使 用 怀柔 太阳 观测 基地 的 Ha 全 日 面 像 进行 测试 ， 测 试 结果 正常 。 此 外 ， 这 个 
系统 不 但 可 以 用 于 Ho 全 日 面 像 的 云 污 染 处 理 ， 也 可 用 在 其 他 全 日 面 像 的 处 理 中 。 

虽然 统一 计算 设备 架构 的 Thrust 函数 库 提 供 了 对 于 任意 长 度数 据 的 排序 函数 ， 但 需要 先 将 数据 放 
入 向 量 容器 ， 且 无 法 实现 数据 集 间 的 并 行 ， 导 致 排序 时 间 过 长 (单个 数据 集 耗 时 0.5 ms， 总 耗 时 
1.7s) 。 而 采用 双 调 排序 不 但 数据 集 内 部 的 排序 效率 高 于 库 中 的 函数 ， 而 且 可 以 实现 数据 集 之 间 并 行 
( 受 共享 显存 的 限制 ， 单 次 可 并 行 数据 集 有 限 ， 经 测试 每 次 并 行 12 个 数据 集 总 耗 时 最 短 ， 为 26. 8 ms) 。 
虽然 该 算法 对 所 排 数 据 的 长 度 有 约束 (长 度 必须 为 2 的 索 ， 因 此 将 极 坐标 角度 设 为 512) ， 但 并 不 影响 
后 续 的 判断 。 

在 云 污染 的 修复 环节 ， 生 成 巴特 沃 斯 低 通 滤波 器 用 时 45.6 ms ， 相 对 于 除 傅 里 叶 正 反 变换 的 其 他 
运算 来 说 ， 这 是 一 个 十 分 耗 时 的 运算 。 主 要 原因 是 每 个 线程 所 承担 的 运算 量 较 大 ， 而 图 形 处 理 器 单个 
处 理 单元 性 能 不 及 中 央 处 理 器 内 核 。 曾 尝试 将 该 滤波 器 做 成 一 张 图 像 保 存在 硬盘 中 ， 让 系统 每 次 从 硬 
盘 中 读 取 然后 复制 到 图 形 处 理 器 中 参与 运算 。 但 测试 发 现 ， 读 取 用 时 77 ms， 再 复制 到 图 形 处 理 需 用 
时 14. 2 ms， 总 耗 时 约 为 在 图 形 处 理 器 中 生成 的 两 倍 。 

然而 ， 系 统 也 存在 一 些 不 足 之 处 。 首 先是 检测 可 修复 云 污染 算法 ， 针 对 云 污染 对 称 均 匀 分 布 的 情 
况 检测 失效 (当然 这 种 情况 极 少 出 现 ) 。 其 次 是 频 域 滤波 的 缺陷 : (1) 当 云 污染 比较 严重 时 ， 如 果 通 过 
增 大 滤波 器 的 截止 频率 更 大 限度 地 去 除 云 污染 ， 则 同时 会 去 除 很 多 太阳 活动 细节 : (2) 图 像 在 太阳 边 
缘 处 存在 较 大 梯度 差 使 得 滤波 时 出 现 一 些 边缘 效应 。 

要 强调 系统 的 主要 目的 是 在 观测 时 实时 显示 云 污 染 修 复 后 的 太阳 图 像 ， 因 而 着 重 关 注 处 理 速 度 。 
虽然 对 修正 后 的 数据 也 进行 了 保存 ， 但 如 果 需 要 进一步 的 对 这 些 云 污染 图 像 进 行 高 精度 处 理 ， 如 对 各 
种 太阳 活动 现象 进行 自动 识别 、 分 割 和 测量 ,采用 文 [5] 提 出 的 多 尺度 形态 学 滤波 更 为 合适 。 如 果 将 
多 尺度 形态 学 滤波 并 行 化 作为 后 续 需 要 考虑 的 工作 ， 同 时 也 需要 考虑 对 由 于 不 精确 的 平 场 改正 、 快 门 
效应 以 及 图 像 饱 和 等 原因 所 致 畸变 的 Ha 全 日 面 图 像 的 修复 ,使 得 系统 有 更 广 的 适用 范围 。 进 一 步 验 
证 系统 涉及 的 云 污染 识别 与 修复 原理 是 否 同样 适用 于 太阳 局 部 像 的 云 污染 处 理 ， 如 果 适 用 再 利用 图 形 
处 理 器 技术 开发 相应 的 处 理 软件 。 

致谢 同时 感谢 由 NOS 主持 的 GONG 项 目 和 国家 天 文 台 怀柔 观测 基地 所 提供 的 Ha 全 日 面 图 像 。 
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A Real-Time Image Processing System for Detecting and 
Removing Cloud Shadows on Ha Full-Disk Solar 


Zhang Nengwei, Yang Yunfei, Li Ranyang, Ji Kaifan 
(Computer Technology Application Key Laboratory of Yunnan Province, Kunming University of Science and Technology, 


Kunming 650500, China, Email; flyingwithcloud@ sina.com ) 


Abstract: Sky clouds lead to Ha full-disk solar observed images covered with shadows, which obscure 
the details of solar features. In this paper, a real-time image restoring system is introduced. For detecting and 
removing cloud shadows, and displaying the restoring image in real-time, Graphic Processing Unit (GPU) 
with Compute Unified Device Architecture (CUDA) parallel programming environment is employed to 
accelerate computing in the system. The parallel computing steps include; (1) identifying heavy cloud covered 
images by calculating the ratio of major and minor axes of a fitted ellipse; (2) recognizing recoverable cloud- 
covered images by calculating the symmetric of limb darkening curve; (3) removing cloud shadows by 
Butterworth low-pass filter in frequency domain. The total processing time is about 0. 7s for one cloud-covered 
image; it meets the need of displaying in real-time. The processing time of each operation taken in GPU is 
measured. The results show that Fourier transform, inverse Fourier transform and Butterworth low-pass filtering 
are the most time-consuming part and together they take up 52.9% of the total time. This paper also evaluate 
the quality of restored image, and the conclusion is that the algorithm can effectively eliminate cloud shadows ; 
meanwhile, it has little effect on the details of solar features. The existing problem and probable improvement 
are also discussed at the end of this paper. 


Key words; Cloud shadows; Parallel computing; Detecting and removing; Butterworth low-pass filtering 


